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In a recent communication to the cond-mat archives, 
Suslov [Q severely criticizes a multitude of numerical re- 
sults obtained by various groups for the critical exponent 
v of the localization length at the disorder-induced metal- 
insulator transition (MIT) in the three-dimensional (3D) 
Anderson model (AM) of localization as "entirely ab- 
surd" and "evident desinformation" . These claims are 
based on the observation that there still is a large dis- 
agreement between analytical, numerical and experimen- 
tal results for the critical exponent g. The author pro- 
poses, based on a "simple procedure to deal with correc- 
tions to scaling" , that the numerical data support v f» 1, 
whereas recent numerical papers find v = 1.58 ± 0.06 

!-!■ 

As we show below, these claims are entirely wrong. 
The proposed scheme does neither yield any improved 
accuracy when compared to the existing finite-size scaling 
(FSS) methods, nor does it give v w 1 when applied to 
high-precision data. 

FSS at the Anderson MIT has a noteworthy his- 
tory, reaching a first peak with the seminal papers of 



Pichard/Sarma 
pecially in Ref. 



and MacKinnon/Kramer MM. Es- 
the groundwork for a reliable, nu- 
merical FSS procedure was laid and scaling curves could 
be constructed that proved the existence of an MIT 
in 3D. In these and later studies based on the same 
analysis technique [Q, the critical exponent v, as esti- 
mated from the divergence of the infinite-size localiza- 
tion and correlation lengths £,(W) at the transition W c , 
i.e., £ oc |1 — W/W c \~ u , is systematically underestimated, 
since the divergent nature at the transition can only be 
poorly captured by FSS of data obtained for small system 
sizes and large errors e in these finite-size data. However, 
as more powerful computers became availably in the last 
decade, one observed a trend towards larger values of 
v k, 1.35 [fn|£3| for e< 1%. 

In 1994, high-precision data (e < 0.2%) showed a hith- 
erto neglected systematic shift of the transition point 
W c with increasing system size. Taking this into ac- 
count phenomenologically, v = 1.54 ± 0.08 was found 
0. A subsequent approach by Slevin/Ohtsuki in- 
corporated these shifts as irrelevant scaling variables and 
further allowed for corrections to scaling due to nonlin- 
earities. With higher-precision data (e « 0.1%), they 
found v = 1.57 ±0.04. Further results for, e.g., the AM 
with anisotropic hopping p3|-|l7j] , the off-diagonal AM 
HQ, the AM in a magnetic field @|l|, confirmed 
this value of v within the error bars (see Fig. [l]). Also, 
v is identical for the MIT as a function of disorder or 
energy 18|l9|| . We emphasize that a properly performed 
Slevin/Ohtsuki scaling (SOS) procedure needs to assume 
various fit functions and that the final estimates are to 
be suitably extracted from many such functional forms 
|Jff8|; bootstrap @-§ or Monte Carlo methods [l§-|lj] 
then need to be employed for a precise estimate of error 
bars. 

We have tested the method proposed by SuslovJ lj] first 
with the transfer-matrix (TM) data of Refs. |§0.[9| 



with e < 0.1%; we find Suslov = 1-75 ± 0.17 for the 
anisotropic and 1.55 ± 0.04 for the random-hopping AM. 
The SOS gives v = 1.61 ± 0.07 g§ and v = 1.54 ± 0.03 
IP^P^I , respectively. Using for a second test energy- 
level-statistics (ELS) data [|l6| with e ~ 1%, we find 
^Suslov = 1.51 ± 0.25, whereas SOS gives 1.45 ± 0.2 
fU-jfl . Last, for artificially generated data with precisely 
known W c = 16.5 and varying v £ [0.5, 2.0] the results of 
the Suslov method are comparable to the results of the 
MacKinnon/Kramer FSS and slightly less reliable than 
the SOS. 




FIG. 1. Results for W c and v, for the anisotropic AM 
|5[[rj using SOS of TM data (open symbols) and ELS data 
(filled symbols) for various fit functions. The error bars show 
the 95% confidence intervals. The accuracy of TM localization 
lengths data is an order of magnitude higher than that of the 
ELS data and the system sizes of TM data are larger than 
for ELS data, giving systematically larger v values for the 
former. The goodness of a fit is reflected in the size of the 
symbol. The 2 thick error bars mark high quality ELS fits for 
large system sizes. The gray o and □ and the corresponding 
error bars (dashed lines) represent t'susiov of TM data and 
ELS data for the anisotropic AM, respectively. The solid line 
marks the result of p). 

We conclude that the method proposed by Suslov also 
yields v « 1.58 and not v » 1 for the MIT of the AM. 

In principle, the Suslov method does not need to as- 
sume any functional form of the FSS curves just as the 
MacKinnon/Kramer method. As a numerical tool, the 
Suslov method is not unreasonable, but certainly not 
better than the established methods: it does not take 
into account the systematic shift due to irrelevant scal- 
ing variables, it relies on an a-priori knowledge of W c and 
inherently produces rather large error bars for the critical 
exponent. We note that Suslov in his numerical test ||] 
used data for 3 small system sizes 6 3 , 12 3 and 28 3 , while 
currently sizes - 50 3 (for ELS) and 18 2 x 10 8 (for TM) 
are standard. It is evident to people with experience in 
FSS that Suslov's erroneously small v is due to his use of 
too few and too small system sizes. 
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In conclusion, high-precision numerics with error e ~ 
0.1% together with all the above mentioned FSS methods 
produce a critical exponent v w 1.58 > 1 for 3D. The 
numerical values of v for dimensions 2 < d < 3 p2~f and 
4 ]22] , |2"c| remain valid, they are certainly not "entirely 
absurd" although there is only limited agreement with 
the field theoretic approach || . Similarly precise data are 
much harder to obtain for our experimental colleagues, 
but recent advances in this direction show a clear trend 
towards increasing v |2^j24]| . 
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